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We examine the Zeno and anti-Zeno eflects in the context of non-Markovian dynamics in entangled 
spin-boson systems in contact with noninteracting reservoirs. We identify enhanced non-Markovian 
signatures in specific two-qubit partitions of a Bell-like initial state, with results showing that the 
intra-qubit Zeno effect or anti-Zeno effect occurs in conjunction with inter-qubit non-Markovian 
dynamics for a range of system parameters. The time domain of effective Zeno or anti-Zeno dy- 
namics is about the same order of magnitude as the non-Markovian time scale of the reservoir 

C'i ' correlation dynamics, and changes in decay rate due to the Zeno mechanism appears coordinated 

with information flow between specific two-qubit partitions. We extend our analysis to examine the 

>— ' , Zeno mechanism-non-Markovianity link using the tripartite states arising from a donor-acceptor-sink 

^^ ' model of photosynthetic biosystems. 

Qh| PACS numbers: 03.65.Yz, 03.65.Xp, 03.65.Ta, 03.67.Mn, 42.50.Lc, 71.35.-y 
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I. INTRODUCTION 



Quantum dynamical semigroups [l| operating via trace preserving, completely positive (CP) maps [2|-y] with inbuilt 

divisibility [5|-l7[ , are well known to provide tractability when characterizing the evolution dynamics of memoryless 

quantum systems. These semigroups constitute the non commutative counterpart of classical Markov theories [Si], 

"S ' and are based on contractions of density operators embedded either in the Hilbert space of Hilbert-Schmidt operators 

Cd ! [l| or the Banach space of trace-class operators [3, [i3| , and involve the forward composition law. Recently, studies of 

^ ' the non-Markovian dynamics of evolution in open quantum systems [6|, |7|, I 111 . Il2l | have attracted increased attention 

O"^] as the Markovian model breaks down in the strong system-environment coupling regime or when un-factorized initial 

conditions exist between the system and environment |13l4l7j , and for which a statistical interpretation of the density 

,-H ' matrix is not guaranteed. In these instances, there is violation of complete positivity of the reduced dynamics of 

^ the quantum system depending on mathematical mappings known as the assignment or extension maps |14l . llSl . Il9| . 

^^ These maps provide description as to how a subsystem is embedded within a larger system, and can assume negative 

CT^ , values for some correlated states. The quantum system may then evolve with signatures of non-Markovian dynamics 

_' ' with novel features of decoherence and non-invariant states even when an equilibrium point is reached. The non- 

i ' invariant states here refer to states which occupy a subspace that is in a state of transience, arising from the presence 

\l of dissipative elements. 

^— ^ Due to the varyin g ch aracteristics of different distance measures (e.g. trace distance, Bures distance, Hilbert- 

• , Schmidt distance (J, |20|), there is no single quantifiable measure, or a unique definition of non-Markovianity in 
. . ' quantum systems. In general, measures of non-Markovianity are based on deviations from the continuous, mem- 
oryless, completely positive semi-group feature of Markovian evolution. Both the concepts of divisibility [3, 0| and 
distinguishability [6|] have been used to define non-Markovianity, and in particular the trace distance Q, D[pi,p2\ 
» ^ is a well accepted metric measure of distinguishability of two quantum states pi , p2, that can be used to check the 

C^ > violation of the complete positivity during the time-evolution of a quantum system. During the time-evolution of a 
quantum state under the trace-preserving CP map p{0) i->- p{t) = $(i,0)/3(0), the trace-distance does not increase 
with time and D[pi{t),p2{t)] < D[pi{0), p2{0)]. Such a monotonically contractive feature is typical of a divisible 
Markovian mapping on the operator space, and the distinguishability between any states is bounded by its value at 
the initial state. An increase of trace distance during a time intervals is thus taken as a signature (sufficient but not 
necessary) of the attributes of non-Markovianity. 

One area in which non-Markovian signatures has potential applications is in the quantum Zeno effect f2ll|, which 
describes the retarded time evolution of a quantum state subjected to frequent measurements. In the limiting case of 
infinitely frequent observations, a scenario that is physically unattainable on the basis of the Heisenberg uncertainty 
[22i | , further decay is inhibited and the time evolution of the state comes to a standstill. The opposite effect which 
leads to enhancement in time evolution is known as anti-Zeno effect [2j]. The theoretical formulations of the Zeno 
effect, which arises naturally from the quadratic short-time behavior of quantum evolution have been advanced 
considerably by Misra and Sudarshan [2J] and Facchi and Pascazio [25l - [27| . The concept that a measurement process 
can decompose the total Hilbert space into partitions of orthogonal quantum Zeno subspaces has been examined via 
the adiabatic theorem [2^ [2a], as an extension of the Misra and Sudarshan theorem [24|. Under this scheme, the 
initial state is confined to single-dimensional invariant subspace, and its survival probability remains unchanged over 



a period of time. Accordingly, different outcomes are eliminated and the system evolves within the same subspaces 
confined by the total Hilbert space. However it is also possible as in Schwinger's non-selective measurements j28| . 
where projections onto multidimensional subspaces are involved, that a quantum system may be steered away from 
its initial state [29||. 

In this paper, we examine the Zeno dynamics of quantum systems from a quantum information perspective, by an- 
alyzing the behavior of dynamical semigroups during occurrence of Zeno and anti-Zeno effects. There are motivations 
to pursue this line of investigation due to the central role played by the decay process during Zeno/anti-Zeno effect. It 
is known that in unstable systems, the occurrence of both effects depend on critical parameters such as measurement 
frequencies and environmental noise [301. Several quantum systems exhibit a combination of Zeno and anti-Zeno 



effects [3l| such as the nanomechanical oscillator [S^ , disordered spin systems [33| and localized atomic systems [3J] . 
The reversal in decay, which is apparent when the system shifts from displaying the anti-Zeno effect to the Zeno 
effect, necessitates some degree of reversibility during quantum evolutions. It is currently not clear as to the precise 
nature of the entities involved during these exchanges, however such reversals may be linked to feedback mechanisms 
within the open quantum system-reservoir model. This has implications for the semigroup law and contractions in 
the Hilbert space of Hilbert-Schmidt operators [l], with likely violations of forward time translations. 

While the evolution within the projected (Zeno) subspace may be unitary, the overall dynamics appears complex due 
to non-Markovian dynamics. The latter feature may be exploited to explore potential applications in quantum control 
[35l ISa ] and information processing 37] involving periodic measurements, which can produce either a Zeno effect or 
an anti-Zeno effect, depending on the parameters of the system. The possibility of protecting quantum states via 
the Zeno effect has been shown in earlier works |38l - l40J . In this work, we highlight the appearance of non-Markovian 
signatures during the Zeno and anti-Zeno dynamics in specific two-qubit partitions of the spin-boson-reservoir system, 
with a broader aim of demonstrating that the context of non-Markovianity can reveal an alternative phenomenological 
perspective of the quantum Zeno effect. 

Due to advancement in femptosecond based spectroscopy, there is currently great interest in the quantum processes 
which underpins quantum coherences and the exceptionally high efficiencies of energy transfer observed in light- 
harvesting systems [4l|-|45| . Accordingly, we examine the occurrence of an environment induced quantum Zeno effect 
associated with dissipative photosynthetic sinks. The sinks form an integral component of the reaction centers (RC) 
in light-harvesting systems [41] by acting as regions where chemical energy is generated. In this study, we redefine the 
role of photosynthetic sinks as detectors of the excitation that is propagating through the system, and examine the 
impact of the Zeno effect on energy transfer mechanisms,. Such effects are not thoroughly explored in biochemical 
systems, and the physical insight obtained from such studies is expected to provide useful information related to the 
quantum control of energy propagation process in artificial light-harvesting systems. 

This paper is organized as follows. In Section [Hi we present several basic elements involved in the dynamics of 
open quantum systems on qualitative terms, with emphasis on the role of assignment maps in determining the nature 
(completely positive, positve and negative) of dynamical maps. A brief description of the pathological elements in 
density matrix operators due to initial system- reservoir interactions is also provided in Section [Hi In Section IIIII we 
present a review of Sudarshan and Misra [24)] formulations of the Zeno dynamics of quantum systems. In Section [IVI 
we introduce a criteria to detect non-Markovianity in Zeno dynamics, and which will be used to analyze information 
flow in the spin-boson-reservoir system (described in Section|V|. Numerical details of the inter-qubit non-Markovianity 
during intra-qubit Zeno/anti-Zeno effect are provided in Section [VTj Using a redefined non-Markovianity measure for 
a tripartite state, the Zeno mechanism in the photosynthetic reaction center (RC) is examined in Sectior fVlII We 
present our main results and conclusions in Section IVIIII 

II. DYNAMICS OF OPEN QUANTUM SYSTEMS AND ASSIGNMENT MAPS 

In the standard approach involving open quantum systems to model dissipative and decoherence processes with 
memoryless time evolution, it is convenient to represent quantum states via normalized density matrices, p, which 
are hermitian (p = p^), and positive definite. The positive definite attribute ensures that density matrices possess 
eigenvalues that are non-negative (0 < Xi < 1) and sum to unity (Tr(p)=:I) in the complex vector space of the Hilbert 
space, 'H. Density operators have the advantage of being able to provide information on the fraction of ensemble 
systems which exist in a given state, when a large number of systems are under study. If p is representative of a 
pure ensemble, then p^ — p, and in the case of a mixed ensemble, we obtain < Tr(p^) < 1. To be identifiable in 
experimental setups, density matrices have to satisfy the completely positivity condition in which positive states are 
mapped into positive states. For states which can be expressed in a factorized density matrix form of a tensor product 
state with environment at initial time i=0, a dynamical evolution to the final state occurs in a tractable form. Here 
the link between initial and final state of the system can be specified by a semigroup of completely positive dynamical 
maps that ensure that the state operators retain a probabilistic interpretation. The possibility of occurrence of a 



collection of (not uniquely) Kraus operators Q , is only guaranteed by completely positive maps so that the associated 
physical process can still continue to be described by CP maps 

p' = Aip) = J2K,pKl (1) 

a 

where the Kraus operators Ka satisfy the completeness relation, ^a-^a^a = 1, and the map is trace preserving. 
Using the properties of a semigroup of completely positive dynamical maps, a quantum Markovian master equation 
of the Lindblad form, and which describes the time evolution of the reduced open system states has been obtained as 

MM ■■ 

^p{t) = -z[H,p(t)]+X^ A,(p(i)) (2) 



dt 



k,l = l 
d 



-i[H,p{t)] + ]^Y.aM (2xkp{t)x\ - {xlxh P{t)}) , (3) 

-i[H,p{t)] + ^7fc Ukp{t)Cl - \{ClCu, p{t)}) , (4) 



fc=i 



where the Lindblad operators, Cki ( Eq. ^) generate the map from the initial to the final density operators of p. 
H arises from a combination of the isolated system Hamiltonian, Jig and a system-environment interaction operator. 
{Xfe}fc=o form the basis in the linear operator space, with xo — ^ denoting the identity. The terms (aki) in Eq. ([3]) 
constitute the positive definite d-dimensional Hermitian Gorini-Kossakowski-Sudarshan matrix A 46] , with a spectrum 
characterized by the decay terms, {"fk}- The first term in Eq. ([3]) (or (j?])) represents reversibility in system dynamics, 
while the symmetrized Lindblad operators are denoted by £k, in which both an operator and its hermitian conjugate 
are labeled by k. The latter incorporate environmental effects within the Born-Markov approximation and therefore 
act as the source of non-unitary dynamics. 

While the Lindblad form in Eq. Q ensures the positivity of density operators at any time, it is stringent in 
being applicable only to weak system-reservoir coupling when Markov approximation holds, greatly simplifying the 
mathematical structure of the mapping procedure. Eq. ^ is therefore not representative of situations when complete 
positivity of density operators is not a necessary feature. There are positive maps where not all eigenvalues of the 
mapped density operators are positive, and positive states continue to be mapped into positive states. And there are 
non-positive maps in which the positive states are mapped into negative states with at least one negative eigenvalue. 
As pointed out in the Introduction, the latter may arise when the system and its environment are initially correlated. 
Pechukas [U] first raised the possibility that non-Markovianity may result as an artifact of the product of the initial 
conditions, Ps{0) €5 Pr{0), where ps {pr) denote the density operator specific to the system, s (reservoir, r). In this 
regard, the Agarwal-Redfield (AR) equation of motion is seen to violate simple positivity in the semigroup form of 
Lindblad for certain initial states j47„ ^] . 

The Agarwal bath [431 ^^^ ^-^^o noted to produce spurious effects in the form of density matrix negativity for 
a range of initial conditions at low temperatures [49[. Such violations may occur in the initial stage of evolution 
dynamics due to backflow of information from the reservoir bath at short times comparable to the bath memory time 
[50j . The Lindblad form fails in these situations and it may be appropriate to use a non-Lindblad set of relations 
to describe the system dynamics during the initial period of quantum evolution, in which the finite time-scale of 
the vibrational environment becomes relevant. In this regard, the non-perturbative hierarchical equations of motion 
(HEOM) technique [5l|, l52| which incorporates higher vibrational energies in the bath as well as a finite time scale of 
the dynamics in the vibrational environment may provide a suitable alternative, as it is dependent on a second-order 
cumulant expansion that is exact for a harmonic bath. 

There has been much debate on the occurrence of the non-positivity attribute in dynamical maps [la. uA. [iSl [l9l [53| . 
with justifications provided |19l l53j to show that maps need not necessarily be in a completely positive form to 
describe the evolution of open quantum systems. The physical interpretation [13, [S^ of non-positive maps has been 
linked to mathematical mappings known as assignment maps |l4| . However there still seem to be a view that any 
physical evolution of reduced density operator MUST preserve positivity, and a violation of the complete positivity is 
considered to be unphysical, for which quantum trajectories lack a physical interpretation. These views are consistent 
with Alicki's [l^ stand that "the complete positivity of the reduced dynamics should be and can be preserved." 
However this reasoning has been refuted by Pechukas [19| on the basis of improper use of "assignment maps" outside 
the weak coupling regime in Ref. [l8| . Non-positivity features are of relevance in experimental techniques involving 
quantum process tomography in systems that possess initial correlations with the environment. The dependence of 



positivity of the map on the interplay between the assignment map and the system-environment couphng has been 
recently shown [54j] . with conditions for positivity noted to arise from correlated system-environment states. Using the 
concept of assignment maps, an earlier work [5^ examined the physical consistency of states for which the dynamical 
map is not positive. Such states were seen to be not amenable to the partial trace operation of extended systems which 
include specific correlations. The term "compatibility domain" was introduced ^53i] to describe the set of states for 
which non-positivity is physically valid. In forthcoming sections, we use numerical results to show the occurrence of 
non-Markovian dynamics in two model systems which incorporate Zeno dynamics. The time domains involved during 
effective Zeno or anti-Zeno dynamics are reasonably matched with the non-Markovian time scale of the reservoir 
correlation dynamics. In the next Section we provide the basics of the Zeno effect. 

III. ZENO DYNAMICS AND NON-MARKOVIANITY 

The Kraus operators in Eq.([T]) may be taken as measurement operators, and generalized measurements may be 
represented by a set of linear maps, where the action of each individual map acting on the density matrix may not 
preserve its trace J55| . Here we adopt Sudarshan and Misra [241] formulations of the Zeno dynamics of quantum 
systems, in which measurement in the Hilbert space H are implemented via a Von Neumann projection operator V 
(which may be one- or higher dimensional) with the Hilbert space range, 'H-p. An initial density matrix po of system 
S in H-p = VH satisfies po = 'Ppo'P, T'^Po'P] = 1- In the absence of any measurement, the state evolves as (24 26] 

pit) - Uit)poU\t) (5) 

U{t) == exp{-in*t) 

where Ji* is a time- independent Hamiltonian with a lower bound. The probability that the system remains within 
H-p is given by P{t) — Tr (U (t) poU^ {t)V) . In the event of measurement at time r, the density matrix p{t) transforms 
as p{t) — -Att VU{t)poU^t)V, where the survival probability within the subspace oi H-p appear as |24l.l26l| 

p{t)=Tt{V{t)poVHt)) (6) 

and V(t) = VU{t)V. In the case of measurements taken at time intervals r = i/iV, the survival probability of the 
quantum state, p{t/N) after N measurements appear as (c.f Eq.®) as 

p{t/N) = Tr {v{t/N)''po[V{t/N)f^) (7) 

At very large N, transitions outside Hp are prohibited, and p{jj) -^ 1, so that the monitored system persists in 
the original state giving rise to the Zeno effect. Misra and Sudarshan [241 showed the semigroup properties of the 
operator V{t) = 'PU{t)P a.t N —> cx), at real time, t. 

The projection-operator partitioning method introduced by Feshbach [57|, divides the total Hilbert space of "H into 
two orthogonal subspaces generated by V and its complementary projection operator 2 = 1 — 7'. We note, following 
Eq.(I7]), that there exists a probability q{t/N) = (1 —p{t/N)) that the quantum system has made a transition outside 
Hp to its ortho-complement, Hq. We consider that the subspace of Hp is spanned by (|l)p,|0) ), while that of 
Hq is spanned by (|1)„,|0) ). Accordingly, |1) denotes the "survived" state due to measurement, while |1) is the 
"un-survived" state. The measurement related decoherence due to leakages between subspaces may be associated 
with quantum jump operators involved in the transfer of states from Hp — >■ Hq. The action of these random jump 
operators form the basis of the quantum trajectory approach [S^, where the density operator is derived from an 
ensemble average of a range of conditioned operators [58| at the select time t. 

We consider an initial state which starts from Hp at t =0, which is applicable to the composite state, \ipt) — 
|1) |0) . After TV measurements during a time duration, i, the system evolves according to the quantum trajectory 
dynamics model [58| 

\A) = VW/N)W,\0), + ^/W/N)\0)J,W,, (8) 



where ^Jp{t/N) {^Jq{t/N)) is the probability amplitude associated with the state residing in the subspace generated 
by P (Q). To simplify the analysis, we assume that the concurrent occupation of states in the two subspaces is not 
viable during the evolution in Eq.®. While a trace preserving map is inherent in the unitary time-evolution between 
measurements in Eq. (jS]), non-Markovian dynamics may arise when leakages between different subspaces are taken 
into account. This possibility is an important consideration in this work, and in the analysis of the quantum Zeno 
effect in the context of non-Markovian evolution dynamics. 



IV. CRITERIA TO DETECT NON-MARKOVIANITY IN ZENO/ANTI-ZENO DYNAMICS 

Under a completely positive map, there occurs no improvement in the distinguishability of a set of states. Therefore 
a completely positive, trace preserving dynamical map of a Markovian evolution on the operator space, p(0) i— >■ 
p{t) = <i>(t, O)p(O), is monotonically contractive with respect to the trace-distance, D[pi,p2] = ^\\pi — P2\\- Here 
||j4|| = Tr[v]4tA] between two states pi, p2. The relation 

D[pi{tf),p2{tf)]<D[pi{U),p2{U)],ioT tf>U (9) 

is useful in identifying violation of the monotonically contractive characteristic of divisible Markovian mapping on the 
operator space. We note that Eq.(j9]) may also imply, with the incorporation of an intermediate time tm-, 

D[pi{tf),p2{tf)] < D[pi{U),p2{t^)], for tf > U,tf > t^ (10) 

tm may not necessarily be the same as U, and Ea. (|10p still applies to a completely positive, trace preserving dynamical 
map. 

To establish a criteria for detecting inter-qubit non-Markovian dynamics, we consider two initial states which are 
subjected to measurements of varying time duration. The trace-distance between these two states is then monitored 
after increasing the measurement duration of both states by the same amount. To simplify the numerical analysis, 
and without loss in generality, we choose two initial states pi(ii) — Pi(0), P2(U) = P2(ji), one of which is subjected to 
measurements of the ideal time duration of r — >■ 0, and is therefore stagnant, while the other initial state is measured 
in intervals of n (which can be taken as variable). For the final states, we select pi{tf) = Pi{t2), P2{tf) — P2{ti -fT2), 
which differ from their respective initial states by being subjected to measurements with interval duration that is 
increased by T2 ■ The specifics of this model is convenient for analysis as will be shown in Section [V] In order to 
determine the occurrence of non-Markovian dynamics, we employ a difference function dependent on trace-distances 

A(n,T2) = i?[pi(r2),p2(ri+r2)]-i?[pi(0),p2(ri)], (11) 

Positive values of A(ri, T2) marks the presence of non-Markovian processes due to underlying information flow from 
the environment back to the system, however we note the different context in which "non-Markovian" is used here 
(reliant on violation of the complete positivity attribute) as compared to a similar term adopted by the chemical 
physical community. In the latter case, the information flow from the environment back to the system, which is 
absent in the Born-Markov quantum master equation (QME), may appear within the Markovian quantum master 
equation by eliminating the Born approximation. 

It is also important to note subtle differences between the terms, "violation of the complete positivity" and "non- 
Markovianity" , as the occurrences of violations of a completely positive map serve only as a sufficient but not necessary 
signature of non-Markovianity. Moreover, non-violation of the positivity may not even indicate the presence of a 
Markovian evolution. This highlights the complex links between dynamical maps and non-Markovian dynamics. We 
next describe the spin-boson system which is used as a framework to demonstrate non-Markovian features of memory 
effects based on Eq. ljlip . and to seek differences between the Zeno and anti-Zeno quantum dynamics. 

V. ZENO DYNAMICS IN THE SPIN-BOSON SYSTEM 

The spin-boson Hamiltonian provides an exemplary dissipative model to examine interaction between a two-state 
subsystem and a thermal reservoir [59| . Here we consider the density matrix of a spin-boson system associated with 
the Liouville equation -r^ = — i[_ffT,p(0]i where the total Hamiltonian H^ = i/qb + Hqs + ^qb-os and _ffqb of the 
two- level qubit assumes the form iJqb = h (^- cr^ + Actx). The Pauli matrices are expressed in terms of the two 
possible states (|0) , |1)), a^ = |0) (1| + |1) (0| and ct^ = |1) (1| - |0) (0|. Afl is the biasing energy while A is the 
tunneling amplitude. 

Each spin qubit is coupled to independent reservoirs of harmonic oscillators, Hos — X^n ^"^q ^q ^q- ^q ^^'^ ^q ^^'^ ^^^ 
respective creation and annihilation operators of the quantum oscillator with wave vector q. The qubit-oscillator inter- 
action Hamiltonian is linear in terms of oscillator creation and annihilation operators i/qb-os = X]q ^ (^q + ^q) '^2- 
The term A denotes the coupling between the qubit and the environment and is characterized by the spectral density 
function, J(a;), which we assume to be of the low- frequency form 

J(^)=5:A^<5(^-^,) = -^, (12) 

q 



The coupling strength between the qubit and the environment, A, is a product of the reorganization energy and a, 
the characteristic frequency of the bath, a is related to the reservoir correlation time, a — j-, and we assume that a 

is small compared to the qubit energy difference, Afl. With the choice of energy unit, Afi^l, -^^ yields a measure 
of a coupling strength that is dimensionless. 

We consider an initial state of the qubit with its corresponding reservoir in the vacuum state at equilibrium, 

|0i) — |l)s <8) rife^i \^k)r — |l)s <8i |0)j. where \1)^ (|0)s) is the upper (lower) level of the qubit, and |0)j. implies that 
all N' wavevector modes of the reservoir are unoccupied in the initial state. During a measurement, a transition of a 
qubit from its excited state to ground state is triggered, which we consider to follow the simple mode of evolution as 
in Eq.®, \4>i) proceeds in time as 

\c^,)^u{t) \l)JO)^. + v{t) |0)J1),, (13) 

Thus we consider an evolution in which the excitation is present in one of the two subspaces as in Eq. (|S|) . In order to 
keep the problem tractable we consider that \1)^. denotes a collective state of the reservoir, |l)r = ^jt^t ^^ A{„}(i)|{n}) 
where {n} denotes an occupation scheme in which there are rii oscillators with wavevector fc = i in the reservoir and 
we define the state \{n}) as \{n}) = |no,ni,ri2...ni..njv'). In the case of ideal measurements, the functions u{t) and 
v{t) in Eq. p^ must always satisfy the relation u{t)'^ + v{t)'^ — 1, however this relation may be violated in the case of 
measurements which introduce significant level of disturbance to the system being monitored [60| . 

The survival probability associated with N measurements (see Eq. ()13p for a single measurement) performed at 
regular intervals r, P^ [t) = u{t)'^^ = cxp(— 7(T)t) where t — Nt. The effective decay rate, 7(r) at small values of r 
is obtained using (A=2, fi=l) Kofnian and Kurizky's formalism [23| 



7(t) = / dujJ{uj)Fr{uj- l:\n), (14) 

Jo 

Facchi et. al. ^6l!l made comparison of the effective decay rate 7(t) to the natural decay rate (without measurement), 
and attributed a smaller 7(t) than the natural decay rate to the Zeno effect while a larger 7(t) (than the natural 
decay rate) is linked with the anti-Zeno effect. The transition point occurs at the intersection of the effective decay 
rate and the natural one at a specific time known as jump time. Ea.(|14p is dependent on the convolution of two main 
functions: (a) the modulating function Fr{io — Ai7) and (b) spectral density J(ij-'). The modulating function is given 
by 



F^(w- Afi) = — sinc^ 
27r 



{u - Ar2)T 



(15) 



and is associated with measurements at intervals of r. In the limit r — > 0, Fr{uj — AO) — > 0, 7(t) — > which explains 
the emergence of the quantum Zeno effect that is independent of the spectral density J(w). 

A third term, f{uj) that incorporates rotating and counter- rotating terms [62] introduces an overall shift in the 
effective decay rate, without greatly affecting further analysis of non-Markovianity with respect to the system pa- 
rameters. We thus employ f{uj)—l (the rotating wave approximation, RWA) to simplify the numerical analysis. We 
obtain the decay rate at the limit, r — >■ cx3, as 70 — j(t — >■ 00) ~ J(Af2). The quantum Zeno (anti-Zeno) effect occurs 
when 7(r) < 70 (7(t) > 70). The parameter range of t and the a, the characteristic frequency of the bath, for which 
these two distinct effects occur is shown in Figure [TJ The quantum Zeno effect always occurs at very short t < 1, 
however beyond r > 1 the characteristic frequency of the bath appears to influence the occurrence of the anti-Zeno 
effect at a < 0.25. 

While the decay rate itself is dependent on the coupling strength term A, the occurrence of Zeno or anti-Zeno effect 
is independent of A, this is partly due to the use of the rotating wave approximation. The influence of A in determining 
Zeno or anti-Zeno dynamics may resurface with removal of the latter approximation, however it is expected to be 
weak compared to the influence of a, the characteristic frequency of the bath. This may be attributed to the non- 
Markovian times-scales at which Zeno dynamics operate. The results here (to be shown in Sectior fVH) appear to 
be influenced strongly by a, which controls the range of time scales to which the bath can respond, and which is 
intricately linked to measurements taken at intervals of r. The results in Figure [T] therefore highlight the critical role 
of a in systems where the RWA can be justified. Moreover these findings suggest the relevance of the Zeno/anti-Zeno 
effect to inter-qubit non-Markovianian dynamics rather than intra-qubit non-Markovianity, in structures with large 
network size and connectivity as is the case with the light-harvesting biochemical systems. We therefore examine the 
link between the inter-qubit Non-markovianity and the intra-qubit Zeno effect in a qubit pair system, next. 




FIG. 1: Occurrence of the quantum Zeno effect or the anti-Zeno effect, based on the ratio of decay rates lili, as a function of 
measurement time duration r and a, the characteristic frequency of the bath. The spectral density function, J(io) in Eq.([T2]) 
is numerically evaluated using a unit system in which h—1, A = 2 and Af2=l. The quantum Zeno effect is marked by regions 
with negative values while the anti-Zeno effect is denoted by positive values. The occurrence of these effects are independent 
of the coupling strength, A (using RWA). 



VI. INTER-QUBIT NON-MARKOVINITY AND THE INTRA-QUBIT ZENO EFFECT IN A QUBIT 

PAIR SYSTEM 



We now extend the composite system in Eq. (jS]) to examine the joint evolution of a pair of two- level qubit systems 
in uncorrelated reservoirs via the following Bell-like initial state 
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where j=l,2 denote the two qubit-reservoir systems associated fmiction Ui{t) in Eq.(|13p. a,b are real coefficients 
and satisfy, a^ + b^ = 1. Using Eq. (jl3|) and tracing out the reservoir states we obtain a time-dependent qubit-qubit 
reduced density matrix in the basis (|0 0) , |0 1) |1 0) |1 1)) which evolves with time duration r as 
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where h = a' + b'v,iT)'v2iT)\ h = abu.ir) u^ir), h = b^v^irfu^irf, h = b^u^irfv^irf, U = b^u,{Tfu2{Tf . 
We assume that the usual unit trace and positivity conditions of the density operator p^^ ^^ are satisfied. The reservoir- 
reservoir reduced density matrix p^^ ^^ is similarly obtained by by tracing out qubit states. Each non-zero matrix term 
of p^i ,2 is easily obtained from the corresponding term p^^^^it) by swapping Ui O Vi. Both matrices possess the 
well-known X-state structure which preserve its form during evolution. Using Eq. (fT3)) and tracing out the one qubit 
and reservoir state, we also obtain the inter-system qubit-reservoir density matrix p ^ ^^ (t) of the same form as in 
Eq.dUl), but one in which /i = a^ + 62yi(T)2«2(T)2, /g = abuiir) V2{t), ^ = &2^,i(T)2«2(r)2, h = b^ui{Tyu2iT)^, 

U = b^U,{TrV2{T)^. 

In Figure [21 we have plotted the trace-distance difference function A(ri,T2) (see Eq. ([TT|) ) as a function of ti and 
T2 (using the qubit-qubit density matrix) at increasing values of a, the characteristic frequency of the bath. We 
recall that ti is inversely proportional to the measurement frequency of one of the initial state, while T2 is the 
increase in measurement time duration of the final states. We note the enhancement of inter-qubit non-Markovianity 
with increase in a, (i.e. increased Zeno effect) at constant coupling strength, A. These results in combination with 
the earlier results obtain in Figure [TJ show that the intra-qubit Zeno effect occurs in conjunction with inter-qubit 
non-Markovian dynamics. A decrease in decay rate is linked with reservoir feedback for the qubit-qubit subsystem. 

In the case of the inter-system qubit-reservoir density matrix p ^ ^^ (i) , we note that the inter-qubit non-Markovian 
dynamics is erased at ri < 0.1 and enhanced at higher ti > 2 with increase in a (see Figure[3]). The results in Figure[3] 
indicate increasing feedback into the qubit-reservoir subsystem, consistent with the intra-qubit anti-Zeno effect that 
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FIG. 2: Sign of Trace-distance difference function Sign[A(ri, r2)] (see Eq.((TT}) as a function of ti and T2 evaluated for the 
reduced density matrix corresponding to the qubit-qubit subsystem (Eq.([T7|)). Sign[x] yields -1, 0, or 1 depending on whether 
X is negative, zero, or positive. Values of a, the characteristic frequency of the bath is indicated at the top of each figure. 
a=b=-j=, coupling strength A = 0.01, number of measurements A^=20. Regions of non-Markovianity are shaded purple, while 
regions which are shaded green denote occurrence of Markovian evolution dynamics. Non-Markovianity increase occurs in 
conjunction with the Zeno effect for the qubit-qubit subsystem. 



applies for the specified range of a shown in Figure [TJ The combined results of Figs [2] and [3] indicate that decrease 
(increase) in decay rate due to Zeno (anti-Zeno) effect may be linked to information fiow between specific two-qubit 
partitions. 

While the results obtained here are specific to the choice of the Bell-like initial state provided in Ea.([T5)). the 
output related to the degree of correlation or anti-correlation of the inter-qubit non-Markovinity with respect to the 
intra-qubit Zeno effect, will vary according to the choice of the initial state. The robustness of the connection between 
non-markovinity and Zeno effect is based on the area size of the non-Markovian region (i.e. the sign rather than the 
magnitude of the Trace-distance) which is evaluated numerically (coloured purple in Figures [U [3]) . The determination 
of the size of this region as a function of the amplitude a, coupling function a, and number of measurements N is 
expected to provide further information on the non-Markovian dynamics, but will be pursued elsewhere. 

We tested the links between the Zeno/anti-Zeno effect and inter-qubit non-Markovianity for all possible two-qubit 
partitions possible within Ea. ([TB|) . and noted a definite dependence on the partition choice (qubit-qubit, reservoir- 
reservoir, qubit-reservoir) in relation to non-Markovian dynamics. We have focussed on the two main subsystems 
which showed significant contrasting modes of information fiow with respect to changes in decay rate arising from 
Zeno or anti-Zeno effect. The results obtained here show that the bath memory time (as seen in changes due to a) 
plays a critical role in the non-Markovian processes of qubit-bath interaction, and may provide a basis for the Zeno 
effect-non-Markovianity link. There is approximate matching in the time domains involved during effective Zeno or 
anti-Zeno dynamics and the non-Markovian time scale of the reservoir correlation dynamics. We bear in mind that 
these results are obtained on the basis of the RWA and the assumption of a constant coupling strength, A. In the 
next Section, we consider a tripartite qubit state constructed from a simple model of the light-harvesting complex, 
accordingly we adopt slightly different criterias for determining the Zcno-effect and non-Markovian regions. 
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FIG. 3: Sign of Trace-distance difference function, Sign[A(ri, T2)] as a function of n and T2 evaluated for the reduced density 
matrix corresponding to the inter-system qubit-reservoir density matrix p^j r2(*)- Values of a, the characteristic frequency 
of the bath are indicated at the top of each figure. a—b=A=, A — 0.01, number of measurements N—20. Regions of non- 
Markovianity are shaded purple, while regions which are shaded green denote occurrence of Markovian evolution dynamics. 
Areas of non-Markovianity increases with a at higher ri > 2, while presence of non-Markovianity is erased for n < 0.1 for 
increasing a. 



VII. ZENO EFFECT AT DISSIPATIVE SINKS IN THE PHOTOSYNTHETIC REACTION CENTER 



It is well known |63l - l65| that several chemical and physics systems can be modeled using a reaction coordinate 
linked to an effective potential energy function with two distinct minima points, and which is amenable to analysis 
in terms of exchange of quantum correlations between the system under study and its immediate environment. The 
light harvesting complex is an important system where the observed long lived coherences, lasting several picoseconds, 
may be linked to quantum information processing features inherent in the propagating exciton (or correlated electron- 
hole pair). The biological pigment- protein complex, called FMO (Fenna-Matthews-Olson) trimer complex in the 
sulphur bacteria [41| constitutes three symmetrically equivalent monomer subunits. Each monomer unit is generally 
modeled as a network of eight bacteriochlorophyll (BChl)a molecules within a matrix of protein molecules. The 
series of excitonic exchanges between molecular sites is described by an exciton Hamiltonian of the form in Eq. (j4]) . 
Local Lindblad terms which represent dephasing and dissipation processes linked to the surrounding environment are 
incorporated during the exchange mechanisms. The exciton energy is irreversibly transferred to a reaction center 
(RC), where it undergoes conversion into chemical energy. 

Earlier works [4J] examined the fine interplay between quantum coherence and enviromental noise (both Markovian 
and non-Mar kovian) to achieve optimal functionality in photosynthetic systems. In a separate study [66|, the non- 
Markovian processes were seen to be prominent in the reorganization energy regime, also noted by Silbey and coworkers 



67| who showed that the reorganization energy and the bath relaxation rate played critical roles during the energy 
transfer process. Chen et. al. [681] presented results which showed the importance of incorporating the site energy 
based electronic coupling correlations. The increased oscillations of entanglement in the non-Markovian regime of 
intra-qubit systems in an earlier work [69] further highlighted the importance of feedback mechanisms in the FMO 
monomer system. However in all these works, there was focus only on the intra-qubit Markovian attribute, and to 
this end, the contributions from inter-qubit Markovian dynamics of photosynthetic systems need further study. 

An important reason to consider inter-qubit Markovian dynamics is associated with the significance of multipartite 
states that act as sources of local and non-local correlations in large photosynthetic membranes which feature in 
many FMO complexes [69|. The intra-qubit Markovian (and non-Markovian) model, which is suitable for the FMO 
monomer, is less applicable in the examination of quantum correlations and exchanges across a wider networked 
molecular system with capacity to hold a large number of entangled excitonic qubits. The latter system may arise 
when entangled photons are absorbed by the photosynthetic membrane, setting up fruitful conditions for a system of 
initially correlated chromophores. In the FMO monomer, it was shown recently j43l] that the unique location of the 
eighth chromophore (at site 8) gives rise to strong inter-monomer interactions, facilitating excitation transfer between 
monomers of the FMO trimer. Hence our approach focussed on inter-qubit non-Markovian dynamics is best suited 
to the realistic condition of an entire photosynthetic membrane constituting many FMO complexes and thousands 
of bacteriochorophylls. In Section IVll we noted the time domains involved during effective Zeno mechanism appears 
matched with the non-Markovian time scale of the reservoir correlation dynamics. This points to the relevance of 
the inter-qubit Markovian dynamics in the context of the Zeno mechanism in applications related to light harvesting 
systems. 

Here we utilize the tripartite states arising from a donor-acceptor-sink model to examine the subtle link between 
the intra-qubit Zeno effect and the inter-qubit non-Markovian dynamics. The Zeno criteria used is dependent on 
the indirect action of the dissipative sink on the non-Markovian dynamics of a specific inter-qubit system. The non- 
Markovianity is evaluated using the trace-distance difference measure of a tripartite state, these criterias therefore 
differ from those used in Section |Vl A prototypical energy transfer mechanism in light-harvesting systems is thus 
modeled via a photosynthetic reaction center (RC) consisting of the donor and acceptor protein pigment complexes 
and a third site acting as the phonon-dissipative sink with a continuous frequency spectrum. The energy from the 
acceptor is considered to dissipate into the phonon reservoir over time. The effective Hamiltonian, Ht, which describes 
of the donor and acceptor coupled to the sink consists of an exciton at the donor site d, coupled to second exciton at 
the acceptor a site, which in turn is interacting with its own source of bosonic reservoir r {h=l) 

Ht = H,+Hr + Hi (18) 

Hs = iOd o% ot + uja (J\ a't + V{ai < + a^ a''_) (19) 

N 

iJ, = ^c^fe6[.&fe, Hi = J2i^k<ylbi + AO, (20) 

k k=l 

where ujd (wq) is the exciton resonance frequency at the donor (acceptor) site , cr!f_(cr''_) denotes the raising (lowering) 
operator of the exciton [70] at the donor site, and a'^{a°-^) denotes the raising (lowering) operator of an exciton at the 

acceptor site. The operator bk (bj^) annihilates (creates) a phonon with frequency Wk in fc-th mode of the reservoir. V 
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is the coupling constant between the acceptor and donor excitons and Lpk is the hnear couphng between the exciton 
at the acceptor site and the phonon with frequency w^ in A;-th mode. 

Using Feshbach's projection-operator method [53], the total Hilbert space of Ht is divided into two orthogonal 
subspaces generated by the projection operators, V and Q where 2 = 1 — 7^. For the photosynthetic reaction center 
(RC), 7^ = cr^ crl + cr° at, and Q = X]fe ^I^fc: such that VQ= QV=0. The density operator associated with the qubit- 
cavity system is given by ps{t) — V pV, p being the density operator of the total system. Using a phenomenological 
master equation that incorporates only energy dissipation from the exciton at the acceptor site (with decay rate Ac) , 
and providing for the presence of a reservoir system with a flat spectrum, the reduced density matrix, ps is obtained 
as 



-Ps = i{psH'^-H'^p,) + X,ntps<Jl, 
H'^ = Ud 4 at + [ua + 8 - i^)alat + V(at a\ + a% at) 



(21) 
(22) 



where 8 arises due to energy renormalization, and the decay rate Ac is influenced by the environmental bath parameters 
such as temperature and distribution profile of vibrational frequencies. Eq. (j2ip is equivalent to the Haken-Strobl 
model at infinite temperature, where pure dephasing is accounted for in terms of a classical, fluctuating field in the 
presence of a sink [67]. By setting the resonant condition oj^ = Wq + (5, and assuming the presence of the initial 
excitation to be in the donor exciton state, the multipartite state of the donor-acceptor-reservoir can be expressed as 



l^t) = |l)d®|0)„®|0) 



(23) 



The presence of an exciton at the donor (acceptor) site is denoted by 11)^; (|l)a)j whilst |0)^ ^^) a) denotes the absence 
of the exciton at the donor (acceptor) site. The reservoir state |0),^=J|j, |0)j, denotes the vacuum state, while \\)^ is 
a collective state consisting of singly excited states of the form, |1)^.= -^ X^fc Vk\^k)r where i?o = \J^Ylik IVfcpj ^-nd 
states which are orthogonal to it |13 |. 

On the basis of the quantum trajectory approach t58[, dissipative processes result in the transfer of states from the 
exciton donor d- acceptor a subspace to the reservoir subspace. Eq. pB]) evolves under the action of the Hamiltonian 
in Eq. ^ as 



1^*) = it |1), |0), |0), + -q, |0), |1), |0), + xt |0), |0), |1), , (24) 

where |^tp (|?7tP) is the probability that the exciton is present at the donor (acceptor) site, and \xt^ — 1 — ICtP ~ V\t^ 
is the probability that the collective reservoir state is excited. For the initial condition, ^o=l, ?7o=0, and we obtain the 
analytical expressions, |^tp = g-Aet/2 |^(;osJ7i+ |^sinJ7i] , and \r]t? = e 

217 = (4^2 _ (^)2)l/2_ 



-Aet/2^ 



sin f2i, where the Rabi frequency 
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FIG. 4: Non-Markovianity in the entangled photosynthetic reaction center (RC) for a state of the form given in Eq. (|25p . 
Trace-distance difference function, Z)(t, t) (see Ea. (|28|l ) is evaluated as a function of t and decay rate, Ac for the density matrix 
corresponding to the exciton-donor(l)-sink(l)-sink(2) subsystem (£0. (127^ 1. a=b=-j= and values of t are indicated at the top 
of each figure. Large damping (Ac) associated with higher rate of detection by the dissipative sink, appear to enhance non- 
Markovian dynamics (positive values of Z)(i, t)). The unit system is based on ft = 1/ =1, with time f obtained as inverse of fJo 
(at Ac = 0) 
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Ea. (p4)) can be extended to examine the joint evolution of a pair of qubit systems (denoted) in uncorrelated reservoirs 
via a initial state of the form 






a\tpt)i \ipo)2 + b\ipt)2 IV'o)i 



= 11; 



dl 



|0)al|O>.l 



IV'o)2 = |l>rf2|0)a2|Ol2 



(25) 
(26) 



where i=l,2 denote the two distinct qubits and a,b are real coefficients (as in Eq. ()16|) '). Prominent non-Markovian 
features are seen (Figure U) in the tripartite qubit associated with the donor exciton(l)-sink(l)-sink(2) partition: 
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(27) 



D[p{t),p{t + T)]-D[piO),p{T) 
D[piO),piT)] 



(28) 



to identify violations of the monotonically contractive characteristic features in the tripartite state of Eq.®. The 
results in Figure 2] show that a large damping or decay rate (Ac) associated with the trapping process at the dissipative 
sink leads to an enhancement of non-Markovian dynamics in the donor exciton(l)-sink(l)-sink(2) partition, Pdi,ri,r2 (t) 
fEa. (l?7l) '). The non-Markovian effects is pronounced at small values of r. The two dissipative sinks which act as 
indirect detectors [7l|, appear to induce a anti-Zeno-like effect facilitating information feedback into the specific 
partition considered here. This result is consistent with those in a recent work v21 where it was shown that repeated 
measurements in disordered systems can induce a quantum anti- Zeno effect which enhance quantum transport under 
certain conditions. It is important to note that the Zeno criteria used here is reliant on the detecting role of the 
dissipative sinks, which differs from that used in Section |V] in which a quantitative relation (Eq. ([T4| ) was employed. 

The observation of increased inter-qubit non-Markovian dynamics due to the dissipative sinks in the protypical 
donor-acceptor-reservoir model may be applied to the larger networked systems of photosynthetic biomolecules. A 
quantitative assessment of the contribution from the anti-Zeno-like action of the photosynthetic sink to the efficiency 
of energy transfer is expected to be numerically intense involving matrix sizes larger than the tripartite state in Eq. (j27p 
Such calculations would need to take into account extrinsic factors such as network size and topological connectivity 
present in the molecular structures of multichromophoric macromolecule(MCMM) systems. 

As discussed in Section |TT1 the density matrix in any physically realistic subspace, in any initially correlated system- 
bath state, appear not to necessitate a statistical interpretation. In general, it is desired to employ exact and non- 
perturbative methods to examine general non-Markovian dynamics of quantum systems. However the quantum master 
equation approach ((Eqs. (ll8l) . (|19p )) and subsequent evaluation of the density matrix of high dimensional systems of 
large biomolecular system present expected insurmountable challenges. In this regard, the results of the set of states 
for which non-positivity appears, should be compared with alternative approaches such as the non-perturbative 
hierarchical equations of motion (HEOM) technique [5l|, |52|, with availability of appropriate computational resources 
in future works. It is hoped that the donor-acceptor-sink model used here, will serve as a starting point to comprehend 
the mechanisms under which non-Markovian and dissipative dynamics appear to coexist, with implications for salient 
properties of chemical and biological systems. 



VIII. DISCUSSION AND CONCLUSIONS 



The observations of long-lasting coherence at ambient temperatures |42| in large biomolecular system have yet to 
be satisfactorily explained due to the sheer complexity of modeling large complex molecules against a background 
of noisy processes. Till now, the exact link between the robust coherence times and the time taken for excitation 
energy to be transferred from the antenna to the reaction center in the FMO trimer complex remains unknown. This 
work is a step forward in a qualitative examination of possible coordination between the observed coherence and the 
environment-assisted transfer mechanisms based on the Zeno-effect. 
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Further insight into understanding of the mechanisms in natural systems will provide innovative ways for robust 
control in open quantum dissipative systems. Photosynthetic complexes have been modeled as quantum channels 
[73| in Ref. [7^ , where noise-assisted energy transport have been examined in terms of transmission of classical and 
quantum information, via characterization by population damping and decoherence parameters. The classical and 
quantum channel capacities of photosynthetic complexes were seen to be enhanced by the presence of noise unlike the 
noiseless system with zero capacity limit. 

The role of a noisy environment which acts as a detection medium with Zeno-like effects, could be further expanded 
for applications related to quantum control and selection of quantum channels with capacities of a quantifiable 
amount. The result in Figure IH even if it is obtained for one specific partition, highlights the possibility of quantum 
control of energy transfer mechanisms through concatenated sequences and pulse interval optimization associated 
with dynamical coupling th eory |75l. l76l| . using artificial light harvesting systems. It is known that open- loop quantum 
optimal control techniques [77J may be utilized to verify quantum coherent transport mechanisms in light-harvesting 
complexes. To this end, Zeno-like effects may be incorporated in optimally shaped laser pulses |771] to initialize a 
photosystem, either in localized or extended mode, to discriminate different routes of transport, and to direct energy 
along critical pathways in light harvesting complexes. This possibility needs detailed study to seek effective quantum 
control for viable applications in energy systems, and especially in the regime of optimal energy transfer. 

In summary, using an entangled system consisting of two spin-boson and and reservoir systems , we have identified 
enhanced non-Markovian signatures which can be associated with the Zeno or anti-Zeno effect depending on the 
subsystem under consideration. The results indicate that changes in decay rate due to the Zeno or anti-Zeno effect 
may be linked to non-Markovian type feedback from the system or reservoir. The time domains involved during 
effective Zeno or anti-Zeno dynamics is of the same order of magnitude as the non-Markovian time scale associated 
with the reservoir correlation dynamics. This study highlights that the context of non-Markovianity can reveal a 
different phenomenological perspective of the quantum Zeno/anti-Zeno effects. Consideration of the Zeno effect in a 
model photosynthetic system highlights enhanced inter-qubit non-Markovianity in the tripartite subsystem involving 
two dissipative sinks. Further investigation is needed to seek a rigorous link between the efficient rate of energy 
transfer in light-harvesting systems and the mechanism by which the monitoring action of the dissipative sinks results 
in notable non-Markovian signatures in specific partitions of multipartite states. The latter states hold much relevance 
in large photosynthetic membranes which constitute many FMO complexes |69| , with capacity to hold a large number 
of entangled excitonic qubits, when conditions for occurrence of initially correlated chromophores become favourable. 
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